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ABSTRACT 



The possible existence of black holes has fascinated scientists at least since Michell and 
Laplace's proposal that a gravitating object could exist from which light could not escape. 

^ ■ In the 20th century, in light of the general theory of relativity, it became apparent that, 

G^ ', were such objects to exist, their structure would be far richer than originally imagined. 

Today, astronomical observations strongly suggest that either black holes, or objects with 

^ ' similar properties, not only exist but may well be abundant in our universe. In light of this, 

black hole research is now not only motivated by the fascinating theoretical properties such 
objects must possess but also as an attempt to better understand the universe around us. We 
review here some selected developments in black hole research, from a review of its early 
history to current topics in black hole physics research. Black holes have been studied at 
all levels; classically, semi-classically, and more recently, as an arena to test predictions of 

/\ ■ candidate theories of quantum gravity. We will review here progress and current research 

j^ . at all these levels as well as discuss some proposed alternatives to black holes. 
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1 Introduction 

This paper presents some developments in black hole research from its very early history to 
modern day. Any manuscript undertaking such a task is bound to be incomplete, the subject 
matter being enormous. What is intended here is a coherent, reasonably self-contained (and 
relatively brief) article capturing essential features in black hole history and research at the 
classical, semi-classical, and quantum level. The goal is to give the interested researcher 
or student an overview of some of the research that has been done and is currently being 
pursued in all these areas within a single manuscript. The style is more of a survey than an 
in-depth study and it is hoped the interested reader will find the references useful for further 
information. Given limited space and time, there are regrettably many, sometimes glaring, 
omissions and entire fascinating areas of research had to be left out. It was therefore decided 
that the bulk of the effort go into reviewing a few selected topics in four dimensional black 
holes within the context of the original general relativity theory of Einstein and Hilbert and 
their natural extensions into the quantum realm. A sincere apology goes to the authors of 
works not included here or accidentally missed. Some of the topics are chosen due to their 
lasting impact in the field as can be seen, for example, by the number of papers appearing 
on the arXiv related to these topics, and are not necessarily new. It is hoped that this type 
of review will give researchers in other areas of gravity a brief overview of the phenomena 
that these recurring topics comprise. 

Black hole research has turned from an obscure, almost ignored area of research to one 
of the most studied segments in gravitational field theory. Today it is common to see more 
than a few black hole papers appear on the pre-print archive on a daily basis. It seems that 
the black hole still has many interesting surprises, from the purely classical, to the purely 
quantum, and everything in between. 

The presentation here is done in a somewhat historical perspective. However, the bulk 
of the results focus on more recent developments as there are a number of excellent books 
and reviews from a purely historical point of view ( [1], [2], [3] and references therein). 

In section 2 we give a brief history of black hole research, which dates at least back 
to the 1780s. In section 3 we present the classical black holes and the fascinating research 
that accompanies them to this day. In section 4 we look at semi-classical research which 
also includes a section on Hawking 's amazing result of black hole radiation. In section 5 
we study results from quantum gravity, primarily the loop approach, which is not to imply 
that other approaches are not fruitful. In starting to write that section of the manuscript it 
seemed that with several differing theories, either no justice could be paid to any of them, 
or else the section had to focus on the research occurring in just one of them. It seems that 
loop quantum gravity is closest to the spirit of the rest of the paper and has also produced a 
number of interesting results. 

Of course, black holes would not be nearly as interesting if not for the fact that there is 
now reason to believe that they may well exist (perhaps even in abundance) in the universe. 
We therefore focus on current astrophysical black hole research in section 6 (much of the 
discussion in this section also applies to the possible detection of primordial black holes [4] 
- [6], although many would have evaporated by the present era.) This is perhaps the fastest 
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changing area in black hole research and it seems that the activity in this field will not be 
dying down any time soon. Finally, in section 7 we also discuss some alternative theories 
of collapse which avoid the formation of a singularity or even the black hole altogether. 

2 Black holes: a short history 

The idea that a gravitating object could exist from which not even light could escape seems 
to date at least back to the work of the Reverend John Michell in 1783 [7]. At this time it 
was already known that light traveled at a finite speed from Roemer's studies of Jupiter's 
moon lo [8]. If light behaved like a particle, with finite speed, why then could it not be 
affected by gravity like other objects? 

In 1796 Pierre Laplace, apparently unaware of Michell's work postulated exactly the 
same thing; that a "dark star" could exist from which no light would escape [9]. Both 
Michell and Laplace calculated that an amount of mass M must be present within a radius 
R = ^^TT^ in order for light not to escape from the object. The circumference corresponding 
to this radius was called the critical circumference. As is well known, this value is in 
(surprising) agreement with the value given by general relativity theory. 

Although it was initially believed that these dark stars could be populous in the universe, 
they were later considered to be at best an academic curiosity, as the size or density such a 
body would possess was considered unphysical. Michell originally calculated that an object 
with a similar average density to that of the sun, would need to be approximately 500 times 
larger in order to stop light. Put another way, an object with the same mass as the sun would 
need to possess a diameter of a mere 20 kilometers. As well, studies in the 1800's indicated 
that light was a wave possessing no mass and therefore it was believed that it would not be 
influenced by gravitational effects. 

The situation changed in 1915 when Einstein and Hilbert formulated the now famous 
field equations of general relativity [10], [11], which in the notation of this manuscript are 
presented as ^ : 

It was not long after the formulation of the field equations that astrophysicist Karl 
Schwarzschild came up with an exact solution which described the gravitational exterior 
of a perfectly spherical star [12]. The Schwarzschild solution is probably the most famous 
non-trivial solution of the field equations and it admits the following well-known line ele- 
ment in the spherical coordinate chart: 



t^^lV c^^d^lV ^/^ ^ ^U ■ (tj 



2GM\ 2 dr^ 



ds^ = - 1 ^— dt'^ + ^;^TTT + r"^ dO"^ + r^ sin^ 6 d(t>'^ . (2) 



c^r 



Two things are immediately apparent in (O: (i) A singularity is present when r = 0. 
This singularity was believed to be of exactly the same nature as the corresponding one in 



In subsequent sections we shall be utilizing geometrized units where G = c = 1. 
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the Newtonian theory and was of little concern, (ii) There exists a singularity in the metric 
when r = ^ , exactly the same radius value for which an object in Newton's theory 
yields an escape velocity equal to c. However, even though it was well known that light was 
affected by gravity in this new theory, it was still believed that stars whose radius was less 
than 2GM ^gj-g {Q ]^Q considered as pathological and not existing in nature. 

This period of complacency did not last particularly long. Many scientists at the time 
began to worry that something seemingly unnatural appeared in the gravitational field the- 
ory and appealed to Schwarzschild's constant density interior solution [13]. For a constant 
density sphere of density po and radius a, Einstein's equations along with the conservation 
law T^'^u = yield the following for the (isotropic) pressure: 



P = Po 



Y^l - 2Mr2/a3 _ ^i _ 2M/a 
3^1 - 2M/a - Vl - 2Mr^/a^ 



(3) 



with r < a. In this expression it was noted that as a approaches |M, the central pressure 
becomes infinite. Thus, it was argued, systems with smaller radius-to-mass ratios could not 
exist in nature. This argument, however, relied on the use of an unphysical matter model 
and perhaps could not be trusted. 

Einstein himself was uncomfortable with the fact that the solution to his field equations 
admitted such bizarre structures [1]. He attempted to disprove the existence of such compact 
objects by studying the circular orbits of massive particles in the Schwarzschild space-time. 
As he made the orbits smaller and smaller, the particle velocities increased until they finally 
reached the speed of light when located at a radius of 1.5 times the critical radius. The 
conclusion was that a spherical object could not exist whose radius was less than this one 
since no particle could exceed the speed of light [14]. This study, though perfectly correct, 
neglected radial motions of the particles which are inevitably present for all particles with 
geodesic orbital radii of less than 1.5 times the critical radius. 

In the same year as Einstein's paper was published (1939), Robert Oppenheimer and 
Hartland Snyder performed an extremely difficult and pioneering calculation. They studied 
the gravitational collapse of a spherically symmetric isotropic dust cloud within full gen- 
eral relativity. Their seminal results were published in an article entitled "On Continued 
Gravitational Contraction" [15]. In this study they demonstrated how the dust cloud could 
collapse and how, when viewed from an external vantage point, the collapse would slow 
down and asymptotically halt as the critical radius was reached. They also showed that, to 
an observer co-moving with the dust, the collapse takes place in finite time. Although dust 
is by no means a realistic matter field, Oppenheimer and Snyder's calculations provided the 
most complete argument of gravitational condensation at the time. 

One of the continuing opponents to black hole formation was J. A. Wheeler. Wheeler 
and many other scientists at the time were quite certain that some physical process must 
intervene during the collapse, and that the scenario played out by Oppenheimer and Sny- 
der's idealized calculation would not occur. One proposal was that the nucleons present in 
a realistic structure would, under extreme conditions, radiate away [1]. Interestingly, this 
idea turns out to be partially correct in the paradigm of black hole evaporation. 
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By the 1960s computers were at the stage where a more reahstic matter model could 
be utilized in gravitational collapse calculations. Such a calculation was carried out by 
Colgate and his collaborators at Livermore laboratory [16]. These calculations, although 
still perfectly spherically symmetric, took into account now well known nuclear physics 
processes inside the matter. These calculations indicated that for a star of mass greater 
than approximately two solar masses, collapse into a black hole was inevitable. Similar 
calculations were carried out in the Soviet Union by Zel'dovich and collaborators yielding 
similar results [1]. Such simulations along with the discovery of new coordinate systems to 
describe black holes [17] [18] aided in easing the scientific community's skepticism about 
black holes. In fact, it is well known that Wheeler became a big believer and is the originator 
of the term "black hole". 

With the possible formation of black holes now accepted by much of the scientific 
community, black hole research saw the birth of the now famous Hawking-Penrose singu- 
larity theorems [19], the no-hair theorem [20], and of course, the Kerr solution to the field 
equations with all its interesting properties and peculiarities, and the laws of black hole 
mechanics. 

The laws of black hole mechanics arose from the analysis due to various talented sci- 
entists and were put in their final form by J. Bardeen, B. Carter and S. Hawking [21]. In a 
nutshell, they can be stated as follows: 

0. The surface gravity is constant on any surface corresponding to a black hole event 
horizon. 

1. If an amount of material of mass SM, angular momentum SJ and charge SQ accretes 
into a black hole, the area of the event horizon responds according to 

^ 6A = n6J + <^6Q-6M, (4) 



8ttG 



where k is the surface gravity of the horizon, $ the electrostatic potential at the 
horizon and Q. the angular velocity of the horizon. 

2. The area. A, of the event horizon cannot decrease. 

3. It is impossible to reduce k, to zero by a finite sequence of operations. 

Although some of these laws have caveats, it was not lost on the scientists of the day their 
amazing resemblance to the laws of thermodynamics. Specifically, the first law of black 
hole mechanics and the first law of thermodynamics would be analogous if one associated 
entropy with the area and temperature with the surface gravity. It was J. Bekenstein who 
took this analogy most seriously and today black hole entropy is commonly associated 
with his name. (For an interesting survey of black hole thermodynamics see [22], [23] and 
references therein.) 

Some research at this time started to focus on quantum properties of black holes. A 
theory of quantum gravity is notoriously difficult to come by, although today there are 
some serious candidate theories. In the late 70s and 80s, in the absence of a full theory of 



Developments in Black Hole Research: Classical, Semi-classical, and Quantum 7 



quantum gravity, scientists started to study black holes from a semi-classical perspective. 
That is, the geometry of space-time was treated classically, but the matter fields propagating 
on the space-time were treated quantum mechanically. The fields were usually "test fields" 
in much the same way as test particles are used. They are quantized on the background 
space-time using techniques similar to many of those employed in quantum field theory 
in Minkowski space-time. There are, however, some issues, ambiguities and subtleties in 
curved space-time that are not present in the corresponding theory in flat space-time. Some 
of this will be discussed in later sections. Excellent expositions of this subject may be found 
in [22] and [24]. 

One very important result that has emerged from semi-classical studies is that black 
holes do indeed, as Wheeler suspected, evaporate. This black hole evaporation was first 
suggested for rotating black holes by Zel'dovich and later, in 1974, Hawking quantitatively 
discovered that all black holes must radiate [25]. This was confirmed and extended by 
D. Page and W. Unruh [26] [27]. To this thermal radiation one could associate a tempera- 
ture, which for the Schwarzschild hole is given by 

^^"^ " 87rGMk ' ^^^ 



as well as an entropy. 



S»»j^A. (6) 



with A being the area of the black hole's horizon. The temperature is very small, approx- 
imately 10^^ K for a black hole of the order of a solar mass. The entropy, on the other 
hand, is very large, being of the order of 10^^ for a solar mass black hole (in J-K^^). As we 
will see, theories of quantum gravity may explain the origin of this large entropy from the 
gravitational degrees of freedom associated with the horizon. 

Since these results, the arena of semi-classical black hole research has been a very 
fruitful one. However, a more fundamental problem remained; a full theory of quantum 
gravity was still (and in many ways still is) elusive. 



3 Classical Black Hole Research 

3.1 Review of important solutions 

Here we briefly review some of the important classical black hole solutions and their prop- 
erties. These solutions are amongst the most studied metrics in general relativity theory. 
A detailed exposition on the mathematical aspects of classical black holes may be found 
in [28] and [29]. 
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3.1.1 The Kottler-Reissner-Nordstrom black hole 

Perhaps the most famous solution to the gravitational field equations is the Schwarzschild 
metric. With charge and cosmological constant it yields the line element: 



/ 2M Q2 A 2\ , 2 dr"^ 2 ,^2 2 ■ 2 /> , 2 



V r r2 3/ 1_2M,Q2_A2 

r r^ 3 

(7) 
Horizons exist where gu = 0. 

It is actually not difficult to derive this solution. Consider a general, static, spherically 
symmetric line element: 

ds^ = -e^(^) dt'^ + e"^'") dr^ + r^ dO^ + r^ sin2(^) V ^ (g^ 

Expression ([8]) yields the following, from the field equations ([T]): 

e""^'') 1 

G*i = — (1 - ra(r),^) - — = SvrT* (9a) 

G; =^^ (1 + n{r),r) - ^ = SttT; (9b) 



-a(r),,7(r),, ) = SvrT^g = SvrT^, (9c) 



For simplicity the cosmological term and the charge term are set to zero, however, it is 
straight-forward to implement these, by including them as part of T^- 
Equation (|9al ) may be utilized to give the following: 

g-a(r-) ^ ^ l\J^2 (rT.ttJ^\ ^j ,-,_.-, 2m(r) 



r 



Ur'f (tI{t')) dr' + 1 =: 1 - ^^ . (10) 



Since the system of equations is under-determined, two functions may be prescribed. Since 
the Schwarzschild solution, if considered in the domain < r < 00, corresponds to the 
gravitational field of a point mass of mass M, we can postulate a stress-energy tensor for a 
point mass with the following Tq component: 

n(r) = -^5(r). (11) 

The r dependence is motivated by dimensionality arguments and neglecting the factor of 
47r would simply correspond to a rescaling of the mass. There is some arbitrariness on what 
the other function to be prescribed can be. However, to satisfy junction conditions implied 
by the above equations supplemented with the conservation law, T\. should be continuous. 
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and therefore set to zero. The remaining unknowns may be solved for by straight-forward 
manipulation of the field equations and conservation law: 



g7(0 =e-«(r)e?«o^ (12a) 

Md{r) 
IGvrr 



K --T^^r . (12b) 



The constant ho can be absorbed into the definition of the time coordinate. 

For the case of a charged black hole (Q ^ 0, A = 0) there are two horizons, one at 
r = r+ = M + \J M"^ — Q^ and another, inner, horizon at r = r_ = M — y^M^ — Q'^. In 
the extremal case, Q = M and the horizons are coincident. 

In closing this sub-section we quote the form of the line element in several other well 
known coordinate systems which historically have shed light on the causal structure (now 
with Q = = A). In Painleve-Gullstrand coordinates, the Schwarzschild metric is cast in 
the form: 



ds"^ = -dP+ldr + J—dp\ + r"^ de + r^ sin^ 9 dif"^ . (13) 

These coordinates are regular at the horizon (although (700 still vanishes there) and readily 
display the no-escape property of this surface. Considering radial light-rays (ds = 0, dO = 
dip = 0), the equation of motion yields: 



dr 




2M 




= ±1- 


\ 


dt 




V r 



Note that for r < 2M the quantity ^ is negative for both solutions, indicating that both 
ingoing and "outgoing" null rays approach r = 0. 

The ingoing Eddington-Finkelstein coordinates are also particularly useful in elucidat- 
ing the "no escape" property of the event horizon. In these coordinates, the Schwarzschild 
metric yields: 

1 \dv'^ + 2dvdr + r'^ d9^ + r^ sin^ 9 dcp^ . (14) 

The nature of these coordinates, along with the causal structure is demonstrated in figured] 
Of course, there also exists the outgoing Eddington-Finkelstein coordinates: 

ds"^ = - (1 - — j du"^ - 2dudr + r"^ d9'^ + r^ sin^ 9 dip"^ , (15) 

which are also illustrated in figure [T] 

Finally we present the line element in Kruskal-Szekeres-Synge coordinates, which 
eliminate the coordinate pathology at r = 2M altogether: 

ds^ = -32 e-''/2A^ dudv + r^ d9'^ + r^ sin^ 9 dip"^ , (16) 

r 
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singularity 




V = constant lines 




w = constant lines 



Figure 1: The ingoing Eddington-Finkelstein coordinates (left) and the outgoing Eddington- 
Finkelstein coordinates (right). The dashed line represents an in-falling time-like particle. 



where r now represents the solution to: 



uv 



1 



2MJ 



\ r/2M 



The Kruskal diagram (also known as the maximally extended Schwarzschild space-time) is 
illustrated in figured The two coordinate patches corresponding to r > 2M and r < 2M 
in metric (|7]) (recall Q = 0, A = in this discussion) are capable of describing regions I 
and II respectively of the maximally extended Schwarzschild space-time. 

Further discussion of horizon-regular coordinate systems may be found in [30]. 



3.1.2 The Kerr-Newman-de Sitter black hole 

This space-time describes a charged rotating ring or, in its asymptotic regime, the space- 
time outside a rotating charged star. It is the rotational analogue of the Kottler-Reissner- 
Nordstrom black hole. Unlike the previous case, here there is no Birkhoff's theorem [31] 
guaranteeing uniqueness of this solution as the exterior of a realistic rotating star. 

The solution, without charge or cosmological constant, was discovered by Roy Kerr 
[32] in 1963, although several others had attempted to find such a solution before him. 
Kerr's original metric yielded the line element: 



ds' 



1 



2Mr\ 



p^ 



I (^du + a sin 6 difY 



+ 2 (du + a^ sin^ 9 dip^ i^dr + a sin^ 9 dip) 
+ p"^ {d9^ + sin'^ 9 dif"^) , 



(17) 



Developments in Black Hole Research: Classical, Semi-classical, and Quantum 1 1 



r=2M 



r=2M 



V t=0 




'' = constant 



t = constant 



-V t=0 



Figure 2: The maximally extended Schwarzschild space-time. 



where 



2 2,22 

p := r + a cos i 



(18) 



with a the angular momentum per unit mass and M the black hole mass. 

With the addition of charge (Q), and in the presence of a cosmological constant (A), the 
line element is expressible in the following formidable form, utilizing the Boyer-Lindquist 
coordinates: 



ds' 



-4;^ (A, - Aea' sin^ 9) dt' + ^ dr^ + ^ dO^ 

P'^2^'^ ' Ar Ag 



+ 



1 



y92S2 

2a 



An (r^ + a^? - Ara^ sm^ e 



sin 6 dip 



[Ae 



r^ + a^ 



Ar\ sin 9 dtdif . 



(19) 



with 



A. 



(r'+a'] 



1 - |rM - 2Mr + Q^ , 
Ae =1 + -Aa^ cos^ 9, S = 1 + -Aa^ . 



(20a) 
(20b) 



This metric is quite complicated to work with and we will therefore focus attention on 
the case (5 = and A = 0. It will also be assumed that a < M. 

In summary, the Q = A = metric possesses the following properties: In the limit 
a ^ this solution goes over to the Schwarzschild solution. The M = limit yields flat 
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Minkowski space-time in oblate spheroidal coordinates. As well, a true singularity exists at 
p = as can be seen from the computation of the Kretschmann scalar 

Rap^sR"'^^^ = -^^8M{r^ - Q? cos^ Q) \p^ - IGr^a^ cos^ Q\ . (21) 

The singularity is located at r = 0, = 7r/2, which in a Cartesian-type coordinate system 
corresponds to x^+y^ = a^, z = 0, indicating a ring-like structure to the singularity. Other 
domains of interest are: 

\)r = r± := M lb \/ Ix-'P — a^; these are the inner and outer event horizons, 
ii) The coordinates t and Lp are not orthogonal. This implies that a geodesic observer "mov- 
ing" in the time direction must necessarily move in the Lp direction. That is, the observer 
is dragged around the black hole in the direction of rotation. This is the famous Lense- 
Thirring effect [33]. 

iii) Related to the previous item, r = r^^ := M it \J M'^ — a^ cos^ d are inner and outer 
ergosurfaces. A time-like observer cannot resist the dragging effects while in this region. 
The light cone structure is studied in, for example, [19], [28] and [34]. 

Figure [3] illustrates the relative locations of these domains, in a Cartesian-like set of 
coordinates. 

There are relatively few treatises on the Kerr geometry, mainly due to the technical 
difficulties involved in working with the metric. The interested reader is referred to [35], 
[36]. 

3.1.3 Exotic Solutions 

Here we briefly describe some solutions which are "exotic" in some sense. By exotic we 
mean solutions that are unlikely to exist in nature but are still of great interest on a the- 
oretical basis. If one allows for a negative cosmological constant, then solutions exist in 
general relativity which can describe black holes with planar, cylindrical, or higher genus 
topology. These are sometimes known as topological black holes. A sufficiently general 
metric to describe such solutions (here with vanishing charge and angular momentum) is 
given by: 

V r J [a^r^ -b- ^^) V ^ 

(22) 
with < if < 2tt. Here, a is related to the cosmological constant via a^ = —A/3, M is 
the mass parameter, and d and b are constants that determine the topology of t, r =constant 
surfaces. The cases are as follows: 

i) 6 = —1, d = —1: In this case constant (i, r) surfaces are spheres (the Kottler solution), 
ii) 6 = 0, lim d = 4: In this case constant (t, p) surfaces are tori. 

iii) 6 = 1, d = 1: In this case constant (t, r) surfaces are surfaces of constant negative 

curvature of genus g > I, depending on the identifications chosen. 

An event horizon exists when (a^r^ — b — ^—) = 0. Such topologies were studied in 
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t — Outer 
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Outer Event 
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Figure 3: Various regions of a Kerr black hole in a Cartesian-type coordinate system. 



detail in [37]. The formation of such black holes from gravitational collapse was studied 
in [38]. 

The metric (|22] ) does not uniquely describe such topological black holes. Lemos and 
Zanchin, for example, have studied a class of black holes which can be cast in the following 
form [39] [40]: 



ds' 



/ 2^2 _ B^\ ^^2 ^2dvdr + r^ {dO^ + V) • (23) 



Here, the values of Bq and the identifications make up either toroidal, cylindrical or planar 
topologies. Specifically: 

i) < < 27r, < (^ < 27r, Bo = V yields the flat torus model. 

ii) — cxD < < oo, {) < Lp < 2tt, Bq = a yields the cylinder, with M the mass per unit 

length and the linear axis coordinate, z, is related to 6 via 9 = az. 

iii) — oo <0<oo, — oo<(/9< oo, i?o = § yields the planar case, with M the mass per 
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unit area. 

This metric differs sUghtly from the metric (l22l) . Gravitational collapse forming such black 
holes was studied in [40]. It can be checked that in all the above cases, the Kretschmann 
scalar blows up at r = 0. Another construction of topological black holes may be found 
in [41]. Research involving the black holes described in this section will be presented below. 

3.2 Developments in classical black hole research 

3.2.1 Quasi-normal mode analysis 

Broadly defined, black hole quasi-normal modes arise from perturbations in some black 
hole space-time. However, the modes are affected by the emission of gravitational radiation, 
which generally has a damping effect on the modes, and therefore these modified modes are 
named quasi-normal ^. The oscillations can in theory be reconstructed via the analysis of 
their corresponding gravitational wave emissions. Reversing the argument, the vibrational 
modes of the black hole are closely linked to the corresponding emitted gravitational wave 
pattern and they are therefore important in light of gravitational wave astronomy. (In the 
case of pure gravitational perturbations, the oscillations are by definition the gravitational 
wave patterns.) It is believed that the gravitational waves due to oscillations produced 
during black hole formation may be strong enough to detect with current or near future 
gravitational wave detectors. The wave signature will be unique, the frequencies tending not 
to depend strongly on the perturbing process and which would yield a direct measurement 
of a black hole's existence and give information on its mass and angular momentum. Two 
or more modes may give useful information such as helping discern if general relativity 
is valid. In particular it may prove or disprove the no-hair theorem of general relativity 
[42]. An excellent expose on quasi-normal modes may be found in the thesis by Cardoso 
[43] as well as the works [44] and [45]. Also, a good reference on non-spherical metric 
perturbations of the Schwarzschild black-hole spacetime may be found in [46]. 

Perturbations of black holes were originally studied by Regge and Wheeler [47]. For 
the Schwarzschild black hole one may, for example, perturb the metric in a way appropriate 
for "axial" perturbations: 



— dt^ + - 
r / 1 



ds^ = -[l — ) dt^ + :; 2p- + r^ dO"^ + r^ sin^ 9 [dip- to dt- q2 dr - qs dOy 



We write 



(24) 
w(r, 9, t) = w(r, e)e'''\ (25) 



and similarly for q2 and q^. 

For perturbations of this form, it can be shown that the system governing the perturba- 
tions can be reduced to a single second-order differential equation, which can be solved by 



■^Here we ignore subtle, but important, mathematical questions regarding the completeness of quasi-normal 



modes. 
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separating the variables r and 6 (for example, see [28]). Briefly, the perturbed field equa- 
tions give a relation between lo and q2 and q^, allowing the elimination of u>. The quantity 
Q := {I — 2M/r)r^ sin^ 6 {d0q2 — drqs) is written as Q = R{r)Q{9) and the resulting 
radial equation is 



r J dr 



1 



2M 



dR 



dr 



fl-2M)i? 



^r 



+ a^R = {), 



(26) 



where /x^ is the eigenvalue of the angular equation and may take on the values /i^ := 
(/ + 2)(/-l)for/ = 2, 3, ... . 

We can make a change of coordinates, 

r, =r + 2Mln(^-l 

so that 



where 



and 



to a Schrodinger-type equation: 




[*? - ^'*^'J 


Zi = -(Ji Zi 


z,M = ^'<'-->. 


V,{t) = 




"/(/ + 1) 6M" 

J.2 ^3 



(27) 



(28) 

The equation (|27] ) is often referred to the Regge-Wheeler equation with (l28l) the Regge- 
Wheeler potential. The above expression is valid for scalar perturbations, the coefficient of 
the last term differing for vector and tensor perturbations (see [48] for details). Later, Zerilli 
derived a similar equation for polar perturbations with a more complicated potential [49] : 



Vi{r)poi 



2(r - 2M) 
r4(nr + 3M) 



[r?{n + l)r^ + SMn^r'^ + 9M'^nr + 9M^] , (29) 



withn := i(/ + 2)(/-l). 

One can compute frequencies for functions which possess the asymptotic form [50] 



Zi{r) 



^Krr, for r=K ^ CXD , 



(30) 



The frequencies are complex and will be denoted as u = cii + ia2- For every value of / 
there exist a tower of modes, denoted here by the integer m. Some of the lowest frequencies 
for the Schwarzschild black hole are summarized in table 1 which hold for both the axial 
and polar perturbations. 

One item of particular interest in modem quasi-normal mode research is that of radiative 
tails. The tail refers to the non-trivial fall-off properties (in time) of the perturbation. This 
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Table 1: A summary of the first four quasi-normal mode frequencies in units of Ma for Z = 2, 3 



M(v 



and 4. The conversion to Hz is given by multiplying the given numbers by 27r-|gi5142Hz. This data 
is from [48] and [51]. 





1 = 2 


/ = 3 


1 = 4 






Mai Ma2 


Mai Ma2 


Mai Ma2 


m = 


= 


0.3737 0.0890 


0.5994 0.0927 


0.8092 0.0942 


m - 


= 1 


0.3467 0.2739 


0.5826 0.2813 


0.7966 0.2844 


m - 


= 2 


0.3011 0.4783 


0.5516 0.4791 


0.7727 0.4799 


m = 


= 3 


0.2515 0.7051 


0.5120 0.6903 


0.7398 0.6839 



phenomenon is sometimes known as black hole ringing, due to the oscillatory behavior of 
the tail. An example of such a tail is given in figure IH for an infalling particle perturbing 
a Schwarzschild black hole. The vertical axis represents the gravitational wave amplitude. 
The dashed line represents an analytical fit using a linear combination of the first two / = 2 
modes. 




Figure 4: A computation of the gravitational wave amplitude of a Schwarzschild black hole per- 
turbed by an infalling particle (solid Une) along with an analytic fit (dashed line) utilizing the first two 
modes corresponding to I ~ 2. Figure courtesy of V. Ferrari and L. Gualtieri, Roma. (From [51].) 
(figure quality reduced for arXiv file size) 



The Kerr black hole is much more difficult to analyze. Some work has been performed 
in [52] - [56]. Studies find that as the angular momentum increases, ai is bounded, whereas 
o"2 is not. However, as a approaches the value of M, a2 tends to zero and the oscillations 
would therefore continue without damping [52]. These modes may not be realized though 
as some studies indicate that the amplitudes of these modes also tend to zero [57]. 

Another interesting phenomenon associated with the Kerr black hole is that of super- 
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radiance. For massless vector and tensor perturbations due to scattering off of the Ken- 
potential, the reflection coefficient can exceed unity if the incoming wave possesses a fre- 
quency below a certain critical value. This is believed to be due to an interplay between 
particle creation in black holes and the Penrose energy extraction process [58] [59]. 

The understanding of quasi-normal modes from black holes is becoming a very im- 
portant issue in black hole physics due to the possibility of detection with gravitational 
wave detectors. There are several astrophysical processes that can give rise to quasi-normal 
modes. For example, an infalling particle could provide such a perturbation. The pertur- 
bation due to infalling extended bodies has also been studied [60]. In this case, the effect 
is smaller than in the point particle case due to interference effects. In the case of fluid 
material orbiting the black hole, modes are only significantly excited when r < AM, and 
is therefore important in the case of unstable orbits. More realistic processes involve large 
scale computations and make up an important area of study in modem black hole research. 
The collapse of a neutron star core has been studied [61] as have the modes caused by thick 
accretion disks [62] and [63]. We will discuss the case of black hole mergers separately 
below. 

3.2.2 Critical beiiavior 

From the point of view of classical black hole physics there are two possible outcomes 
that result from the evolution of regular Cauchy data: Either a black hole forms or it does 
not. In this paradigm, an interesting question to ask is what happens in the regime that 
straddles this bifurcation? This question was originally studied by Choptuik utilizing a 
spherically symmetric massless scalar field minimally coupled to gravity [64]. This choice 
of matter field is convenient as one does not need to worry about possible formation of field 
condensates (stars) in this model. Complete field dispersion or black hole formation are the 
only possible outcomes. Samples of the two outcomes are displayed in figure |5] 




'ft 



M^ 




Figure 5: Evolution of scalar field data undergoing gravitational collapse. In the left diagram a 
black hole forms at late time whereas in the right diagram the field disperses. Figures courtesy of 
M. W. Choptuik, Yukawa International Seminar, (figure quality reduced for arXiv file size) 



The problem was tackled as follows: The initial data contained one tunable parameter, 
usually denoted as p. By changing this parameter the numerical evolution would either 
form a black hole or not. Ideally, the problem is set up so that there is a simple relationship 
between the size of this parameter's numerical value and the "strength" of the gravitational 
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interaction as there is no way a priori to know if the evolution will result in an event horizon. 
Therefore many evolutions need to be run, in each trial changing the value of the parameter, 
before the critical point is reached. 

For the minimally coupled massless scalar, (j), the field equations yield: 

Rfiiy - -j^RQiiv = Svr f (t),^,4>,u - 2<^,a'A'" 9iJ.u j , (31) 

whereas the conservation law gives rise to the wave equation: 

ct>;;^ = . (32) 

With a metric admitting the line element as in dSJl, except that now e'^'^' — > e^^'''' *) =: 6^ 
and e"'^'^' — > e"*^'''*-' =: a^, the field equations go over to: 



Or a? — 1 
a 2r 



"^(h? + ^ 



--0 , (33a) 



^.^ ".^ a^ - 1 „ .--, , 

-i ■ =0 , (33b) 

b a r 

^ - Airrcf) rT =0 ■ (33c) 

' 

The wave equation, after some mild manipulation yields the condition: 

l4>,) =lfr^V) • (34) 

b / ,t r'^ \ a y ,. 

Dirichlet and Neumann boundary conditions were used for the metric and the scalar field. 
At r = 0, the condition of space-time regularity was imposed. As well, the initial (t = 0) 
scalar field profile is also fully specified along with an initial velocity. 

Choptuik found some unexpected and very interesting behavior when studying the 
above system of equations. In the cases that were barely super-critical (i.e. cases where the 
parameter was tuned to values just above the black hole formation limit), Choptuik found 
an unexpected relationship. He found that the masses of the "mini" black holes formed 
followed a scaling law of the form 

M = C{p- pcT ■ (35) 

Here, pc is the critical value of the adjustable parameter. That is the value (within numer- 
ical precision, which was 10^^^ for the original study) p attains when it is exactly at the 
bifurcation point. C is a constant that depends on the particular initial data profile chosen 
and, perhaps most interestingly, the exponent t] possesses a value that is universal for all 
initial scalar field data^ . For the minimally coupled, massless scalar field, Choptuik found 



^This exponent is often labeled as 7 in the literature. However, in this manuscript 7 plays multiple roles and 
therefore rj is used to avoid confusion. 
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that rj PS 0.374 . No matter what the initial data profile, the masses of small black holes 
followed the law (1351 ) with this numerical value of the exponent. 

Since Choptuik's original work, various authors have extended the domain of study to 
other coordinate systems [65], different matter models ( [66] - [75]) and different numbers 
of dimensions ( [73], [76] - [79]). Also, work in non-spherical symmetry has been done as 
well as other couplings and some analytic analysis of the problem ( [80] - [87]). Some of 
the possibilities and results are summarized in table 2. Further, the collapse can be classified 
as type-I, where there is a minimum finite black hole size, or as type-II, where it is possible 
to form arbitrarily small mass black holes. The phenomenon discussed here applies to the 
type-II case. 

Table 2: A summary of some fields studied and their critical exponents. 

Field 7] 

Massless real scalar field 0.37 

Massless complex scalar field 0.39 

Massless complex scalar field (angular momentum) 0. 1 1 

Massive real scalar field (small mass) 0.37 

Radiation 0.36 

SU(2) 0.20 

Gravitational waves 0.36 



The behavior discussed in this section is reminiscent of certain phase transitions in 
statistical mechanics. One example is the liquid-gas phase transition. Near the critical 
temperature, Tc, a substance on its boiling curve will possess a discontinuity in the density 
of the two phases of the form 

pi-Pg^ (T, - T)" , (36) 

where pg and pi are the densities in the gas and liquid phases respectively. Ferromagnetic 
materials obey a similar law near the Curie temperature, T^, 

m^{Tc- T)" , (37) 

where m represents the magnitude of the magnetization vector. 

There is another interesting phenomenon associated with the threshold of black hole 
formation, namely that of self-similarity of the space-time. As this phenomenon is mainly 
of interest on the side of the transition where black holes do not occur (i.e. slightly sub- 
critical) it is beyond the scope of this review and is omitted due to article size considerations. 
As well, we have only scratched the surface in citing the large number of works in this field. 
The interested reader is referred to the thorough reviews in [88] and [89] and references 
therein. 
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In closing this section we quote the words of C. Gundlach [88] on this topic "Critical 
phenomena are arguably the most important contribution from numerical relativity to new 
knowledge in general relativity to date." 

3.2.3 Black hole mergers 

The two-body problem in general relativity is extremely difficult to analyze, mainly due 
to the non-linearity of the field equations. Any, even somewhat realistic, system needs 
to be evolved numerically in the strong-field regime and even with modem computational 
technology it is still a taxing problem. Black holes, being "simple" objects, and having 
the possibility of producing a strong enough gravitational wave signal to detect are natural 
objects to consider in the two body problem. This system provides an excellent arena to 
study the strong-field effects in G. R. without the complications introduced by material 
(i.e. non-gravitational) effects. Since the two black hole system is unstable, it is expected 
that at late times the solution will approach that of a Kerr black hole. The case of a black 
hole-neutron star merger in full general relativity has been recently analyzed in [90]. 

Binary black hole systems were first numerically modeled in 1964 by Hahn and 
Lindquist [91] on a 51x151 mesh utilizing axial symmetry in a head-on collision scenario. 
Within 4 hours the evolution had proceeded 50 time steps in their simulation. A decade later 
Smarr [92] and Eppley [93] also constructed simulations of head-on collisions in the hope 
of studying the gravitational wave emissions. With the announcement in 1990 of the LIGO 
gravitational wave observatories to be constructed, gravitational wave simulations and the 
two body problem were taken up in force. In such studies the numerical methods almost 
invariably involve some form of 3 + 1 split or null variants of it. 

The black hole merger is usually separated into 3 regimes: The inspiral stage, the merger 
stage and, finally, the ringdown stage (see [94] for full details, which includes an earlier 
fourth stage, the Newtonian phase). 

In the inspiral stage, gravitational wave emission has the greatest effects on the dynam- 
ics. It is expected that in the case of small to medium mass ratios, large eccentricities in the 
orbit will decay, yielding a roughly circular orbit by the end of this phase [95], [96]. Ex- 
treme mass ratios are expected to possess high eccentricities [97], [98]. This phase is often 
well modeled by post-Newtonian and higher-order methods. If the mass ratio is extreme, 
the smaller partner may be viewed as a test particle in the background space-time of the 
larger partner. Methods exist to calculate gravitational wave emission in this "test-particle" 
case (see, for example [99] and references therein). 

An interesting point is that the orbit of the small companion in an extreme mass-ratio 
scenario will not generally lie in a plane, due to frame-dragging effects. Therefore, the 
geometry of the surrounding space-time can be well probed by the small companion and 
this information would be transmitted by the gravitational wave emission which is believed 
to lie within the future LISA detector's bandwidth. 

The merger stage is extremely complicated and requires full numerical investigations. 
This is a very short lasting phase, perhaps lasting two cycles of the emitted gravitational 
waves. The luminosity of gravitational waves emitted at this stage may approach the order 
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of 10^^ J/s and the frequency of the gravitational wave approaches that of the dominating 
quasi-normal mode of the resulting black hole. An enormous amount of energy is liberated 
in this phase which totals in the neighborhood of three percent of the rest mass energy of 
the system [94]. 

Finally, the ringdown phase refers to the settling down of the single black hole formed 
to a Kerr black hole. Quasi-normal mode analysis is quite important in this stage as one has 
a perturbed black hole space-time radiating away energy via gravitational wave emission. 
The ringdown is dominated by the following frequencies [42], [94], [100]: 
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(38b) 



with J the angular momentum of the final resulting black hole. The equivalent of one or 
two percent of the rest mass energy is radiated during ringdown. At late time the emission 
is dominated by the radiative power-law tail. 

As discussed above, much work in this area is numerical, due to the complications 
presented by the Einstein field equations. However, there exist constraints in the system 
of equations which must hold throughout any evolution scheme for it to be valid and this 
complicates the numerical evolution. Mathematically, these constraints arise as follows: 
Consider a metric with elements g^j^y and define a unit normal, n^ , to a class of time slices. 
We also define g^y := g^y — n^riy. On a space-like hypersurface the following (Gauss- 
Codazzi) relations hold: 

i?tp. =~9\R\m ~C g% ~9\ - K^Ky, + K^.Kyp , (39a) 

D^Ky„ - DyK^„ = - nxR\p^ T. 9^ ~9l ■ (39b) 

Here R^upa is the curvature tensor constructed with g, D is built with the connection of 
g, and K^y is the extrinsic curvature of the hyper-surface. It should be noted that these 
relations hold regardless of the field equations. By pulling back (I39al ) and (I39bl ) to the 
hyper-surface, and utilizing the field equations, one obtains the constraints for the Cauchy 
problem in general relativity. Therefore, the constraint equations enforce solutions to be the 
allowable data sub-sets permitted within general relativity theory. 

It is a major focus of modern numerical research to find a scheme which ensures that 
the constraints are not seriously violated in the subsequent evolution. As the continuum 
ADM form of the field equations are weakly hyperbolic, this is no trivial task and gener- 
ally some clever method needs to be devised such as consistent modifications to the equa- 
tions of motion. An evolution program often "crashes" or becomes "unstable" when con- 
straint violating modes grow to a size that is considered unacceptably larger than discretiza- 
tion/truncation errors. A breakthrough in the long-evolution stability problem occurred re- 
cently in 2005-06 with the advent of two methods that yield particularly stable evolutions (at 
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least relatively). These are the generahzed harmonic coordinates with constraint damping 
(GHCCD) method [101] [102] and an improved variant of the Baumgarte-Shapiro-Shibata- 
Nakamura method with moving punctures (BSSN) [103], [104], [105], [106]. Briefly, in 
the GHCCD method one adds to the field equations a function of the constraints which 
is designed to minimize or dampen the constraint violation. A recent proposal for such a 
counter-term in the case of the harmonic coordinates was devised in [107]. In the BSSN 
scenario one defines a conformal metric, 

5ii:=e-%, e^ = ~g^'\ (40) 

with g the spatial three-metric in the ADM decomposition, as well as a conformal trace-free 
extrinsic curvature quantity 



ij ■ — ^ 



Kij - -QijK 



(41) 
with K the trace of the extrinsic curvature. The conformal connection is given by 

r:=g''r^k- (42) 

In the BSSN scheme $, Aij, f * and the lapse and shift are considered the basic variables 
of the evolution. The reason this scheme is utilized is that the long ranged degrees of 
freedom can easily be isolated from the non-radiative ones. Another reason is that the 
constraints can easily be substituted into some of the evolution equations thus implementing 
some of the constraints at the dynamic level. Also, with appropriate implementation of 
gauge, the evolution equations are hyperbolic, which is related to the fact that the connection 
(|42l) is treated as an independent quantity. The "moving punctures" refer to the fact that the 
black hole singularities are represented by punctures which move within the grid although 
there is no evolution at the puncture itself. 

An excellent review of the progess in black hole mergers may be found in Pretorius' 
review [94]. We summarize here some recent results. 

For the scenario involving two equal mass black holes with minimal spin and eccentric- 
ities the amount of energy released in the last stages before a Kerr black hole remnant is 
approximately 3.5% of the system's total energy. The resulting Kerr black hole possesses 
a spin parameter of a ss 0.69 ( [94], [105], [106], [108] and references therein). After the 
"collision" the gravitational waveform is dominated by the fundamental of the quadrupole 
moment of the quasi-normal mode of the final black hole. Also, when the flux is near its 
maximum, and subsequently, the waveform may be closely represented as a sum of quasi- 
normal modes, which is surprising as this is expected to be a highly non-linear regime. 

If one removes the restriction of equal mass, but maintains minimal eccentricity and 
spin there is a decrease in the total energy emitted as well as the final spin of the resulting 
black hole. Also, although the quadrupole is still dominant, higher modes become non- 
negligible due to the reduction in the symmetry of the problem. This symmetry reduction 
is also responsible for an asymmetric gravitational radiation beaming. This delivers a recoil 
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to the produced black hole in the orbital plane as there is net momentum carried away by 
the radiation. The maximum velocity imparted due to this effect seems to be approximately 
175 Km/s when the mass ratio is 3:1 [94]. 

For the case of equal mass, nominal eccentricity but non-minimal spin, a new degree of 
freedom is introduced in this case as the individual black hole spins can be aligned in various 
directions compared to the orbital angular momentum. If the net spin angular momentum 
has a component parallel to the orbital angular momentum then more energy will be emitted 
than in the corresponding zero-spin case. If there is a component anti-parallel, less energy is 
emitted. Some particular studies indicate that a pair of holes, each with a w 0.76, radiated 
approximately 7% of their rest mass energy when the spins were aligned with the orbital 
angular momentum. In the case where the spins were anti-aligned, only approximately 2% 
of the rest mass energy was radiated. The final black hole spins were approximately 0.89 
and 0.44 respectively [109], [110]. 

The next case presented is for two equal mass black holes with minimal spin but large 
eccentricity. This is the case studied in the first complete merger simulation by Pretorius 
in [101]. In this study two localized scalar field profiles were initially employed which 
collapsed to form black holes. The outcome of the collapse essentially yields a two black 
hole vacuum. The two field profiles were given equal magnitude but opposite direction 
boosts, with zero boost corresponding to a head-on collision. The final result, merger or 
separation, depends on the value of the single parameter, k, measuring the strength of the 
boost. An interesting result was found. Near the threshold value of the boost parameter, for 
a given class of initial profiles, the number of orbits, n scale as 

e"oc \k-kc\~^ , (43) 

where kc is the threshold value of this parameter. The exponent /3 was found to possess a 
value of approximately 0.34. As it was noted in [94], this behavior is similar to that of test 
particles in equatorial orbit around a Kerr black hole. Those test particles which are near 
the capture threshold approach unstable circular orbits of the Kerr space-time and possess 
a scaling behavior similar to (l43l) . 

Recently, there have been studies providing simple formulas for the calculation of the 
final spin in a binary merger [111], [112]. 



4 Semi- Classical Black Hole Research 

4.1 Review of semi-classical theory 

Before discussing semi-classical relativity, we shall need a result from flat space-time for 
future use. Let us begin by studying the real massless scalar field in Minkowski space-time. 
Such a field obeys the wave equation: 



.M 



'>" - 0. (44) 
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The field can be decomposed into Fourier components of positive and negative frequency 



oo 



(45) 



As is well known, tlie number of particle excitations associated with a particular state may 
be deduced from the number operator: 

N = Y,aian- (46) 

n 

Under proper Lorentz transformations we have (dropping subscripts on w) 



oo oo 

Y, [a„/„(x)e--* + a1,/:(x)e-*] = J] [a:„/;,(x')e— '*' + at;/-(x')e-'*' 



, (47) 



n=0 m=0 



where the prime denotes quantities calculated in the boosted frame. The relation between 
the primed and unprimed modes is 

Jm ^ / _, V^mnjn + PmnJni ! 
n 

Jn = / ^ \pinmJm ~ PnmJm\ i (4°) 

m 

where a and (5 are the Bogoliubov coefficients of the transformation. 
The unprimed and primed vacua are defined via: 

a„,|0)=0, a'^|0'> = 0. (49) 

In particular, the unprimed and primed creation and annihilation operators are related as 



/ ^ O/.nm.O'm ' Pnm^- 



't 
nm"'m 



(50) 



under Lorentz transformations, as can be seen by inserting (1481 ) into (|47] ). Note that if the 
(i coefficients are not zero, |0') / |0) and the two vacua will «ot be the same. That is, 
one observer's zero particle state will not be a zero particle state for a Lorentz transformed 
observer. Another way to view this is that the number operator in the primed frame does not 
agree with the number operator in the unprimed frame. In Minkowski space-time, proper 
Lorentz transformations respect the condition /3„m = and therefore uniformly boosted 
observers will agree on particle content. This is not true for the case of observers which 
are accelerating, even in flat space-time. This leads to an interesting effect known as Unruh 
radiation. Discussion of this is beyond the scope of this manuscript. 

The semi-classical approach is based on the premise that the gravitational field remains 
classical, but the matter content is quantized. The most straight-forward way to incorporate 
this into Einstein's theory is to write Einstein's equations as 

R^.v - \R9t.. = Svr (V| T^Mi) IV") =: Svr (T^,) . (51) 
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That is, the expectation value of a stress-energy tensor operator that depends on quantized 
fields, 4>i, is utilized as a source term in the gravitational field equations. Of course, one 
may still add a purely classical piece to the right-hand-side of the equations, T^,^ , so that 
(^1 T^y{4>i) \ip) yields quantum corrections to the classical theory. The simplicity of equa- 
tion (ISTI ) is deceptive. The matter fields themselves are, of course, quantized on the curved 
space-time, generally leading to a complicated dependence on g^y on the right-hand-side. 
There is also the issue of what state the fields should be in. As well, the quantity {T^y) will 
diverge and therefore a regularization of the effective action leading to {T^y) needs to be 
performed. 

There are several states of particular importance in semi-classical black hole physics. 
These include, but are certainly not limited to, the Hartle-Hawking vacuum, the Unruh 
vacuum, and the Boulware vacuum. The Hartle-Hawking vacuum corresponds to a black 
hole in thermal equilibrium with a bath of thermal radiation, the Unruh vacuum corresponds 
to a state with a particle flux at future infinity, and the Boulware vacuum to a state where 
particles do not traverse to infinity, as measured from near the black hole. 

In the case of a scalar field, the regularization leads to a renormalized value of the 
cosmological constant and a renormalized value of the Newtonian gravitational constant. 
As well, the divergence in the effective action also leads to a term possessing fourth-order 
derivatives of the metric and terms quadratic in the curvature as the divergences possess the 
form: 

(V)d» « I [^9,. + BG,y] + ((i)C(i)i?^^ + i2)C(2)H^^) \na . (52) 

Here {i)H^y and i2)H^y are given by: 

{i)H^iu ■.=2R-fj,-i, — 2gf^yR'^p + -g^yR — 2RR^y , (53a) 

{2)Hfj,u '■='^R,T-u;a ~ ^fj.u'^^p ~ ^^fiuR'fp - 2R^Rau + -g^yR^ Rafi , (53b) 

and a is Synge's world function [113], which is equal to one-half the square of the geodesic 
distance between two nearby points. The limit o" ^ needs to be taken and the logarithmic 
divergence may be removed by renormalization of the constants (i)C and {2)C- Although 
these terms are due to the metric dependence of (T^i/) itself, they are often written on the 
left-hand-side of the field equations: 

Gf,u(ren) + \ren)9fMi^ + {l)C(re„) {l)H^,u + (2)C(^e„) {2)H^,y = 8lT ^T^y + (?),,,) ^^^^^ J , 

(54) 
where the subscript (ren) denotes renormalized values. Equations (l54l ) are sometimes 
known as the semi-classical Einstein field equations. 

An approximation often employed in semi-classical research is a perturbative approach. 
That is, the metric is written as: 

9^lu = {0)9^^ + (^h^u , (55) 



26 A. DeBenedictis 



where tQ\g^u is a classical "background" metric and the small constant parameter e vanishes 
in the limit h —^ so that h^^ represents first-order quantum corrections to the classical 
metric. One customarily limits calculations to order linear in e. The Einstein equations then 
yield an Einstein tensor of the form: 

Gf_iu = {o)Gf_iu + e {h)Gf_iu , (56) 

with (ri)G^v being due to the quantum part of the metric and {q)G^i, satisfies the usual 
Einstein equation for the classical stress-energy tensor. One then constructs a stress-energy 
tensor on the classical background by some means, denoted here as e (q) {T^u), and attempts 
to solve ffijG^u = Stt (q) {T^u) for h^iy. This is the back-reaction problem. 

4.2 Developments in semi-classical black hole research 

Perhaps one of the most amazing predictions to come out of semi-classical theory is the 
evaporation of black holes. The calculation was first carried out in Hawking's paper "Par- 
ticle creation by black holes" in 1975 [25]. Although this is now an old result, given its 
importance it is appropriate to review it in this section. 

We will study this effect in the Schwarzschild black hole ^ in the forms given by (fT4l) 
and (IT5] ). We will also be considering Hawking's original argument where he considered 
a collapsing star with asymptotically flat regions. The Penrose diagram for the collapsing 
star is shown in figure |6] Any direct interaction of the quantum particles with the material 
making up the stellar body will be ignored since the gravitational field of the star is the main 
ingredient for this effect and not the star itself. 



Surface of star 




Figure 6: Penrose diagram depicting the surface of a collapsing star. 



For simplicity a null scalar field will be used to illustrate this effect. At Z^ , there are 
no particles present and the state will be denoted as |0j„). It will be found that, at X+, this 
state corresponds to a state with a non-zero particle content. This is due to the presence of 
the non-trivial gravitational field. 
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Since X^ corresponds to approximately Minkowski space-time, we can define modes 
there with frequency w' and therefore the solutions to (l44l) in Schwarzschild space-time 
may be written as 

'/^L',. = /:'(x)e-*-'\ (57) 

i.e. the standard particle states of flat space-time. The field is written as before, with 
appropriate creation and annihilation operators satisfying a^ |0m) = 0: 

</> = / [aL,/:,(x)e--'- + dlr:,{^)e^'^\ dJ . (58) 

(Note that in this argument we only have inward propagating waves so only inward 
components are necessary here.) 

At X+ the space-time is also Minkowski. We can define modes there as well, which in 
terms of the null coordinate u have the form: 

^u^.u = /(x)^e-'"" , (59) 

which again correspond to the usual flat space-time modes with positive frequency. These 
are outgoing modes. 

Now, we can express ^p in terms of ip' utilizing a similar transformation as in (1481 ): 



^= I \oil,^ip - Puj'^^'*] dv . (60) 

As in the flat space-time case, should the second Bogoliubov coefficient not vanish, there 
will be a disagreement in particle content between an "in" observer and an "out" observer. 
Let us consider the transformation between the in and out observers. For a single out- 
going mode we have: 



e 



-iuj{v-2r*) ^ ^-iuj[v-2r-iM\n{^-l)]_ ^^y-^ 



Consider now tracing back the path of a particle (outgoing) as it skims the horizon, near r = 
2M. If the star is collapsing, the level surface representing the stellar boundary (assumed to 
be near r = 2M) changes with time as the particle traverses through the star. The equation 
describing this surface is given by r ss 2M + kq{v — vq) + 0{v — vq)"^. Adopting the 
eikonal approximation, and assuming most of the change in the eikonal takes place near 
this surface we get 



UJU ~ OJ 



v-4M + 2Ko{v-vo) + ...-^Mln( ^!^i^2_Zl + ... 



(62) 



with V < vq. 

Concentrating on the dominant part of the eikonal, we can see that the ingoing eikonal 
(which we denote as Q) corresponding to the outward eikonal (|62] ) is given by: 

Q,{v) « -AMlo In {{vo-v))+ constant , (63) 
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the constant shift being irrelevant. Comparing this expression with (1601 ) it may be seen 
that, as a Fourier transform, Lpuj,u has both e"*"^ '" and e^*"^ ^ components and therefore the 
con^esponding (3^i^ does not vanish in the transformation. There are therefore out particles 
even in the absence of in particles. The conformal diagram corresponding to the collapse 
into, and subsequent evaporation of, a black hole is shown in figure |7] A much more 
detailed calculation, which includes a derivation of the actual particle spectrum, may be 
found in [114]. 



Surface of star 




Figure 7 : Penrose diagram of stellar collapse taking into account subsequent black hole evaporation. 



This evaporation issue is related to what is sometimes known as the information loss 
problem where "information" falling into a purely classical black hole is presumably lost 
due to the fact that a black hole is described solely by its mass, charge and angular momen- 
tum and therefore the evolution is not unitary. It is thought that if the resulting black hole 
radiation spectrum is not truly thermal then the information may re-emerge from the black 
hole in this way. Later it will be shown that a full theory of quantum gravity may resolve 
this in another way. 

Since the Hawking result, much work has been done in semi-classical gravity, major 
advancements having been made in the 1980s. Much of the early works concentrated on 
methods to obtain sensible, convergent quantities within semi-classical theories. The inter- 
ested reader is referred to [24] and references therein. In the following we limit our survey 
to some issues which directly involve the back-reaction problem as the literature and topics 
in the field of semi-classical gravity are almost as vast as in the classical theory making 
anything resembling a complete coverage nearly impossible. 

More recently, semi-classical black hole research has focused on modeling the pertur- 
bations on the classical background geometry due to quantum fields and their fluctuations. 
Study of the latter effect makes up the arena of stochastic gravity. These problems require 
computation of the renormalized stress-energy tensor on the background, which is then to 
be used as a source for the metric perturbations. 



Developments in Black Hole Research: Classical, Semi-classical, and Quantum 29 



In the black hole context, Hiscock, Larson and Anderson have calculated the back- 
reaction effects of scalar, spinor and vector fields inside a Schwarzschild black hole's event 
horizon [115]. To construct {T^,^)^ they utilized various approximation schemes developed 
in the literature [116] - [120] and the DeWitt-Schwinger expansion [121]. They studied 
the quantum field's back reaction on the anisotropy of the interior as well as the first-order 
correction to the Kretschmann scalar. In summary they found the following: Spinors and 
minimal and conformally coupled scalars tended to decrease the anisotropy as the singular- 
ity was approached whereas vector fields tended to increase the anisotropy. Regarding the 
effect on the Kretschmann scalar, it was found that the minimally coupled massive scalar 
field and spinor fields tend to slow down the rate of increase of curvature as one approaches 
the singularity whereas other couplings and fields tended to increase the rate of curvature 
growth. 

A similar problem was studied for the case of cylindrical black holes (or black strings) 
in [122] utilizing the stress-tensor approximation in [116] for quantum scalar fields. It 
was found there that the conformally coupled scalar also tended to increase the growth of 
curvature near the horizon. In this case, utilizing the cylindrical version of the metric (l23l) 
as a background space-time, with cylindrical (rather than one null coordinate) it was found 
that 



6K 



10 vr TT \M J \ a 



(64) 



where 6K is the perturbation on the background Kretschmann scalar. Here T is the interior 
time coordinate (corresponding to the radial coordinate outside the black string) and the 

horizon is located at T 



{4My/3 



Of course, the above results are only valid insofar as the perturbation is valid and the 
results cannot be extrapolated right down to the singularity. 

In spherical black holes an enormous amount of work has been done in calculating 
field expectation values and stress-energy tensors of various fields and couplings. There 
has also been much work in trying to produce the A/ 4 entropy of black holes from semi- 
classical theory. The amount of work in this fascinating area is too large to even begin 
reviewing here. We refer the interested reader to [159] and references therein along with 
the books [24] [22]. 

For cylindrical black holes the amount of work is much more modest. Very briefly, in 
the context of the cylindrical black holes, Piedra and de Oca have studied the quantization 
of massive scalar and spinor fields over static black string backgrounds [123] [124]. They 
have calulated {T^,^) up to second order in the inverse mass value. Dias and Lemos have 
studied magnetic strings in anti-de Sitter general relativity [125] and the scalar expectation 
value, (</>^), has been computed in [126]. 

Interestingly, critical behavor has also been studied in the context of semi-classical grav- 
ity. For technical reasons, much of this work has been done in two dimensional tensor- 



''From now on in this section we drop the (0) subscript in front of (T^i/) and it is understood that (T^^) 
refers to the first-order quantum correction to the stress-energy tensor and is constructed with the zeroth-order 
(classical) background metric. 
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dilaton gravity as {T^y) may be determined from the trace anomaly along with the mild 
assumption of conservation. Ayal and Piran, for example, obtained a critical scaling expo- 
nent oirj Ki 0.409 [127]. A slightly different model, utilizing a conformally coupled scalar, 
was analyzed in [128] and a critical exponent of 0.5 was found. In [129] yet another model 
was employed which allows one to turn the quantum effects on and off. A critical scaling 
exponent of 77 « 0.53 was found in this study. In [130] the authors calculated the quantum 
stress tensor on a classical background spacetime with perfect fluid source. The quantum 
effects then were treated as perturbations of the classical fluid gravitating system. They 
found that a mass gap exists when r/ > 0.5 so that there is a minimum size to the black 
holes formed. The case r/ < 0.5 could not be studied as the semi-classical approximation 
breaks down in that regime. 

In stochastic gravity one takes the level of quantum approximation one step further, con- 
sidering the effects of the field fluctuations. Given limited space we cannot cover stochastic 
gravity here. The interested reader is referred to the reviews [131] - [134] and references 
therein. 

5 Quantum Black Hole Research 

5.1 Quantum gravity 

Attempts to quantize gravity go almost as far back as the dawn of quantum mechanics. One 
of the earliest arguments for the quantization of gravity, in fact almost as old as general 
relativity itself, is that if T^y is inherently quantum, then so it should be with the gravita- 
tional field which it produces [135]. However it can be argued that this is not necessarily 
the case [136]. Although it is now generally believed that the gravitational field must be 
quantized in some way, there is still some debate on this necessity. 

In the early days of quantum theory, the first person to realize that there would be 
serious problems applying those techniques to gravity seems to have been Matevi Bronstein. 
Bronstein had the insight to deduce that quantum theory could not be applied in any obvious 
way to a theory that was background independent. It was possible, according to Bronstein, 
that the ordinary notions of space and time would have to be abandoned [137]. Other 
pioneers in the field included P. Bergmann and P. Dirac. 

After quantum field theoretic techniques were sufficiently developed, an obvious ap- 
proach to quantizing gravity was implemented as a simple background expansion: 

Qfiv = Vfiv + eh^t, + 0{e^) . (65) 

As is now well known, treating the perturbations as fields on the background metric (rj^^) 
yields a non-renormalizable quantum field theory with divergences commencing at the one- 
loop level for gravity with matter couplings. 

In the early 60s, Feynman, working at tree-level, computed transition amplitudes and 
demonstrated that reasonable results are obtained. This gave hope for this line of quan- 
tum gravity research. However, he noted that at loop level problems began to arise which 
required the introduction of Faddeev-Popov ghosts by DeWitt [139]- [142]. 
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In the 70s it became generally accepted that gravity coupled to matter will be non- 
renormalizable. It was, however, found that one could add a spin | particle to general 
relativity which yields a theory finite at two-loops. Thus began the field of study known as 
supergravity [138]. 

Finally, it was shown explicitly in 1986, by Goroff and Sagnoti that, at two-loop order, 
finite S matrix elements could be attained if the gravity action contained a counter-term of 
theform^ [143]: 



£, 



209 



ct 



{d - 4) 2880(167r)2 



-9R''',sR'',.R'\p 



(66) 



with d the effective regularized dimension of the space-time, thus explicitly showing the 
divergent properties at two-loop level. 

On the canonical side, DeWitt in 1967 publishes what was originally thought of as the 
"Einstein-Schrodinger equation" also known as the Wheeler-DeWitt equation [144]. Some 
argued at the time that the problem of quantizing the gravitational field had been solved. 

This Hamiltonian approach begins with the familiar ADM decomposition of space-time, 
as illustrated in figure [H In the figure qij is the metric of the three-surface S and N and A^-' 
are the usual lapse function, and shift vector associated with the ADM decomposition. 




ds^= -(N--N'NJdt\2N'dXjdl +q,,dx'dx' 

Figure 8: ADM decomposition of space-time into "space" and "time" 



With this decomposition, the gravitational action may be written as: 

1 



16.j''j/'' 



n'^'qab 



+2NMJ'^ (g-i/2n"^) + iVgi/2 (r(^) - q-^U,M"' + l^^^^) 



(67) 



with n := W^Qab, n''^ := g-^/^ [K"'' - K q"''], K := K^^Qab- Kab is the extrinsic 
curvature tensor and over-dots denote differentiation with respect to the slicing time, t. In 



At one loop order the divergence possesses terms proportional to R and ij^^i?'"^ and is therefore finite 
in the absence of matter. 
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this scheme, the field variable is q and its conjugate momentum is 11" . The equations of 
motion go over to 



7V2 



2vf) 


"^-i/2n^/ 




= =: -V^ , 


i?(3) - q-'liab^'-' + ^g-^tf 


= 0=:-S. 


Hamiltonian density: 


Hg = T^n^^g,, - 
iDvr 


Cg = Ni,V 


'' + N S , 



and the following symplectic structure: 

{n'^^(x), g,,(y)} = W7r6l6',^6i^, y), 
{n'^''(x), n^'^(y)} = = {g,,(x), qUy)} , 
which after quantization leads to the famous Wheeler-DeWitt equation: 



1 

—QabQcd + -^Qacflbd 



6 



+ q'/^R^'^ 



^(g) = 0. 



(68) 



Sqac 5qhd 

The Wheeler-DeWitt formulation suffers from some problems. The configuration field (3- 
metric) does not appear as a gauge field. As well, there are inconsistencies with certain 
transition probabilities in the path-integral version. 

There are several other candidate theories of quantum gravity. These include the sum- 
over-Euclidean geometries developed by Hawking [145] and its Lorentzian counterpart, 
the causal set approach of Sorkin [146] [147], dynamical triangulations [148], and other 
theories, including loop quantum gravity. Interestingly, the causal set approach predicts the 
existence of a small positive cosmological constant of the order of that required to provide 
the observed acceleration of the universe. 



5.1.1 The loop quantum gravity program in brief 

Loop quantum gravity provides a promising quantization scheme for general relativity. 
There is a Hamiltonian approach and a covariant approach, yielding a spin-foam model, 
so named due to the resemblance of the Feynman diagram analogs to a foam of bubbles. 
We will concentrate here on the Hamiltonian approach, which is perhaps a bit more per- 
spicuous. For a nice review of the spin-foam approach, the reader is referred to [149]. 

It was noted by Ashtekar [150] that general relativity can be very neatly reformulated 
in terms of a densitized triad, E^- instead of the metric: 



qq 



ElE^d'^ 



with 



E" 



-e''^''eijke\e\, where 



^ a^b"ij- 
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This is sometimes known as the phase-space representation and the indices i, j etc. denote 
the orthonormal components. In terms of the new variables, the ADM action (|67] ) may be 
written as 



/ 



dt 



1 

1 Ti^ . T?h 



E1K\ - NbV'' -NS- N'cijkE^^K 



d^x 



(69) 



KabE j6^^ . In this scheme, the canonically conjugate variables are 



with K'^ := ^_-_ 

E"^ and K\. 

The symmetry group in S is SO (3). We can write K\ in terms of the fiducial so(3) 
connection^ on E, F*^: 

jK^, = A\ - T\ , 

where 7 is known as the Immirzi parameter. The "modified" connection A^^ can be defined 
by this equation. In terms of this new connection the action may be written as 



, dt 
Svr./ ./s 



E'^iA^ - NbVb -NS- N'Gi 



d^x , 



(70) 



I 

Gi := daE'\ + e^j^A^^EI , 
with symplectic structure: 

{E-{^), A\{y)} =87rj5%6)6{^,y) , 
{E-i^), E)iy)] = {^^,(x), A\iy)} = . 

The above formulae do not distinguish between SO{3) and SU{2) and both these groups 
possess the same algebra so it is customary to work in SU{2), the indices i, j now coupling 
quantities to the su{2) algebra. The quantization of this system yields the canonical version 
of loop quantum gravity. £'" and A^^ are the Ashtekar - Barbero variables. The quantum 
versions of the equations of motion yield the quantum Einstein equations: 



Gi\^)=DaE^^\^)=Q, 

E^ab - (1 - ^')K\Gi 
1 



Va\^ 



sm 



Vd^t{E) 






m = 



(71a) 
(71b) 

(71c) 



with 



^ab '--daA f, — ObA a + e ji^A ^A ^ , 

DaVi ■■=daVi - eijkAlv'' . 



*The fiducial connection is that yielded by the solution of Cartan's structural equation, 9[ae\i + 



e^jt.r-'r e u = 0, where e'a is the standard (non-densitized) triad, qab — e^a^^Sij. 
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Note that now there is a constraint equation for Gi (the Gauss constraint). This reflects 
the fact the triad possesses a rotational freedom; one can choose different frames locally 
by rotating the triad. This redundancy is eliminated in the new Gauss constraint. To the 
above system one could add matter couplings by supplementing the action with a matter 
term and quantizing appropriately. The problems associated with (1681 ) are not present in 
this representation. 

Before continuing, we shall make a few comments about this quantization scheme: 

1. The scheme is background independent and respects diffeomorphism invariance. The 
choice of time slicing is arbitrary and does not affect the physics. 

2. A superpartner can be accommodated and therefore supersymmetry can be incorpo- 
rated. This has been done [151]. 

3. Instead of the Einstein-Hilbert action, one can accommodate geometric actions made 
up of arbitrary curvature invariants. The scheme is generally similar to that outlined 
above. 

4. Higher dimensions can be accommodated. 

It should be noted that 2, 3, and 4 are not required but simply can be accommodated. 

What is of interest is the holonomy of A as it is transported around what aie known as 
spin-networks, (first introduced by Penrose [152] and utilized early in loop quantum gravity 
by Jacobson and Smolin [153]) and the state vectors |^) which are functions of this holon- 
omy. The concept of time evolution is now encoded in terms of how the interrelationship 
of the network, which describes space, evolves. The details are beyond the scope of this 
manuscript but the interested reader may find them in [154], [155], [156]. What is of par- 
ticular interest in the context of black hole research is that this theory predicts that on the 
small scale, space is discrete^? Classically, the area may be constructed out of the triad via 



A{S)= / ^Jn-EiubE^d's, (72) 

with the n vectors denoting normals to the 2-surface S. To go over to the quantum theory 
one replaces the classical triad with the corresponding quantum operator. The picture that 
arises is that each fiber of the spin-network that pierces a surface S endows it with a certain 
amount of area and geometry, via the introduction of an angular defect (see figure |9l). 

The triad possesses the following spectrum when acting on the state functions of loop 
quantum gravity: 

E.iSiWiSi) $^' {hM]) = (8vr/27)2 [jp{j^ + 1)] c^i (he[A]) , 

where <!>■' (/ie[^]) are state functions which depend on the holonomy of A, h(.[A], along an 
edge e of the spin-network. The jp are half-integers and Ip is the Planck length. Therefore, 



'it should be emphasized that this is a prediction of the theory and is not put in "by hand". 
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Figure 9: A gravitational spin-network endowing a surface S with area and geometry. The spin- 
network punctures the surface S at p. The surface is flat everywhere except at the punctures. There, 
the spin-network can be pictured as "tugging" on the surface, endowing it with geometry and in- 
troducing a local angular defect on the surface. The n's are the nodes, which are associated with 
volumes. 



for the area operator^ : 

Areas |^) 



SttII^YI 



jpijp + l)\^) . 



(73) 



Notice that the eigenvalues of area are discretel 

One can also construct the classical volume utilizing the triad: 



V 




eabcEfE^.Ele^^'^ 



d^x 



(74) 



This can be replaced by its quantum analog and volume eigenvalues can be calculated. The 
results are somewhat complicated and we omit them here. However, the volume is also 
discrete. 

We will see below that these two operators are of extreme importance in loop quantum 
gravity black hole research. 



5.2 Developments in loop quantum black hole research 

There are a number of results regarding black holes in loop quantum gravity. We shall con- 
centrate here on what are arguably the two most significant results; namely the source of 
black hole entropy and the resolution of the singularity problem. These are of importance 



^We are making an assumption here regarding how the spin-network pierces the surface S. The general 
case yields eigenvalues which are slightly more complicated than ( 173b . 
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because it has long been believed that any viable theory of quantum gravity should ex- 
plain where the enormous entropy of a black hole comes from and it should also eliminate 
singularities present in the classical theory. 

5.2.1 Black hole entropy 

The subject of black hole entropy has been one of intense interest ever since Bekenstein's 
calculations [157]. Many methods have since been utilized to calculate the entropy (see 
[158], [159] and references therein for excellent reviews of the subject). One belief is that 
the source of this entropy is strictly gravitational in origin. That is, one should be able 
to define microstates in a full quantum theory of gravity which, when counted, yields the 
correct entropy law. This has been done within the framework of loop quantum gravity. 

The basic idea is as follows: The gravitational spin-network pierces the surface corre- 
sponding to the horizon of the black hole. As described above, this endows the surface with 
area and geometry (see figure [TOl). The entropy is given by the logarithm of the number of 




Figure 10: A gravitational spin-network giving a spherical black hole horizon its geometry and 
area. 

loop quantum gravity states that give the surface a fixed area, ao- This counting is non-trivial 
as for a black hole of reasonable size there could be an enormous number of punctures, with 
various values of jp. The total area is given by summing up all the contributions from all of 
the punctures. The total area is therefore given by^ 



p 
which is obviously the sum of eigenvalues associated with the edges puncturing the surface. 



'We are making a similar assumption here as in the previous footnote. 
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Notice that this formulation yields the entropy as a function of 7. The quantity 7 therefore 
has to be set by other means or else it can be set by demanding that the leading order term in 
the entropy calculation agrees with the Hawking-Bekenstein result of ao/4. Several studies 
have set the Immirzi parameter in this way [160] - [163]. 

There is another subtlety which complicates the expression which is topological in ori- 
gin. To illustrate this we begin by noting that the action giving rise to the surface states (and 
thus the surface Hilbert space) is a Chern-Simons action: 



7- _ ^'^O rr 
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(76) 



where the trace of the SU{2) connection, A^^ is taken over the su{2) indices. The form 
of this surface term was first calculated in the pioneering work of [164]. Later this was 
generalized to arbitrary genus g surfaces, as above [165]. This term arises from the fact 
that at the inner boundary (the isolated horizon [166]), the triad and the connection cannot 
be fixed independently and are actually related (hence only the connection appears in (1761)). 
The isolated horizon boundary conditions reduce the degrees of freedom and the above 
action can be written in terms of a C/(l) connection [164]: 
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where W represents U{1) connections on the boundary surface, which are restricted by the 
value of the bulk SU{2) connection penetrating the surface at that particular point on the 
horizon. Finally one can construct the symplectic structure on the boundary as was done 
in [164]: 
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(78) 



with k := 



QQ 



^^, _-^. (an integer) known as the Chem-level and with SW and 5'W tangent 
vectors in the space of U{1) connections defined on the horizon. Note that we now have 
a topological U{1) theory on the boundary. The number degrees of freedom are related to 
the number of topologically independent closed paths one can construct on the punctured 
surface. 

In the case of a surface with spherical topology S"^, one can place punctures on the 
sphere and then define a closed path around each puncture. This would seemingly yield 
a 2N dimensional phase-space as each closed path (cycle) also has a conjugate open path 
associated with it (chain) (roughly speaking each cycle represents a configuration variable 
that must have a conjugate momentum, represented by chains). However, these are not all 
independent degrees of freedom. On a sphere, note that going around the A^-th puncture is 
the same as going around all the other punctures but in the opposite direction. If the cycles 
are denoted as rji, this may be expressed as: 



7?1 -772 • ■■■■VN-l =Vn ' 



(79) 



Therefore, the topology reduces the number of independent degrees of freedom. 
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When one quantizes the system, quantum states Vm are obtained for integers m = 
{mi, .., rriN-i) with nii S {1, .., fc} [167]. Note that in this sense, the integers mi play a 
role similar to the magnetic quantum number in ordinary quantum mechanics. The condi- 
tion (1791 ) gives rise to a constraint: 

mi + ... + ruM-i = —mN- ■ (80) 

This restriction is the quantum analogue of the Gauss-Bonnet theorem for a sphere. Note 
that one now has A^ generators and one constraint. Thus, for a spherical horizon, states can 
be labeled with m = {mi, .., mjsf) subject to constraint dSOl ). In other words, not all states 
that yield classical area uq are allowable. Only the ones meeting the condition (l80l ) are to 
be counted. The details of the counting may be found in [168] and references therein. Only 
the results will be cited here. 

A very careful numerical counting of acceptable states for a spherical horizon was per- 
formed by Corichi et al in [162]. Those authors performed the counting with the projection 
constraint (l80l l and without considering it. They found the following: Without considering 
the projection constraint they found that the entropy obeys the ao/4 law, provided the Im- 
mirzi parameter is set equal to 7 ss 0.274. This is the same value found in [161] and [168]. 
With the projection constraint it was found that the number of acceptable states is reduced 
in such a way that does not involve 7. The result is 

S=^-^ln{ao). (81) 

For the case of a genus g surface, the situation is slightly more subtle as the paths around 
the punctures can also be related to paths around the genus holes of the surface. This yields 
a quantum Gauss-Bonnet theorem of the following type [165]: 

Vg+l ■ V9+2 ■ ••• • rig+N = ryiTl'^fSr^ • ••• • VglgVg^lg^ , (§2) 

where r/i through rjg denote the paths associated with the genus holes and rjg^i to rjg^N 
denote with the paths associated with the spin-network punctures. Utilizing this topological 
condition, the entropy of a (j( > 1 horizon is given by [165] 

S=^ + {g-l) [In(ao) - ln(47r7(5 - 1))] , (83) 

provided that the Immirzi parameter is set to the same value as in the spherical case. There- 
fore, the same value of the Immirzi parameter yields the first-order ao/4 term for all cases 
whereas the sub-leading term depends on topolgy and is independent of 7. This behavior is 
consistent with other, non-quantum gravity approaches to calculating black hole entropy of 
g > horizons [169] [170] [171] [172]. 

An ambiguity exists for g = I due to the decoupling of the triad and the connection 
at the horizon in this case. However, one may analytically extend (1831 ) to g = 1 yielding 
an ao/4 entropy without logarithmic correction. This result is consistent with studies of 
g = I horizon entropy utilizing non loop quantum gravity techniques [170], [172]. The 
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g > I result however, is qualitatively different from the g = case, and therefore cannot 
be extended to reliably encompass the g = horizons. This is due to the non-trivial inter- 
play between the spin-network punctures and the genus holes of the surface (note that the 
coefficient of the logarithmic correction in (1831 ) differs from (ISTl) by a factor of 2 for g = 0). 



5.2.2 Removal of the classical singularity 

Another problem that a quantum theory of gravitation is expected to resolve is that of the 
singularities that exist in the classical theory. This is also related to the problem of in- 
formation loss associated with black holes. This singularity issue has been studied in the 
framework of loop quantum gravity in the case of mini-superspace models. That is, models 
where the full system is first reduced to a mechanical system which consists of only the rel- 
evant degrees of freedom. One then quantizes this reduced system. This is done due to the 
technical difficulty involved when trying to work with the full theory. There are currently 
several black hole studies available within the symmetry reduced models [173] - [175] and 
here we shall be outlining the approach in [175]. A related method was studied in [176] 
where it was shown that a Nariai universe replaces the classical singularity. 

The basic ideas are as follows: Construct an evolution equation utilizing the Hamilto- 
nian constraint and check if the evolution remains finite at the point corresponding to the 
classical singularity. Also, one may compute operators and their expectation values which 
encode the information about curvature and which classically diverge at the classical sin- 
gularity. If they remain finite in the quantization the singularity is avoided in the quantum 
theory. 

We will focus on the most studied black hole in quantum gravity, the Schwarzschild 
black hole (A = = Q), whose line element for r < 2M can be written as: 

ds^ = -^^— + ("^ - l) dR" + r2 d02 + t2 gj^2 Q ^^2 ^ (g4) 

where T is the interior time coordinate (corresponding to r in dTJl) and R is an interior 
spatial coordinate (corresponding to t in (|7])). 

Recall that the conjugate variables in loop quantum gravity are the densitized triad and 
the modified SU{2) connection. A pair is is constructed which respects the symmetry of 
the space, R x SO{3) [175], [177]: 



A\Ti dx"" =CT3 dR + (ari + hT2)d6 + (ar2 — 6ri) sin 9 dip + t^ cos 9 dip 

Eyda =PcT3 sin 9dR + (pan + PbT2) sin 9 de + {paT2 - PbTl)d^ , 



(85a) 
(85b) 



where the tj denote the standard su{2) basis. The quantities a, b, c and pa, Pb, Pc are to be 
determined and act as conjugate "position-momentum" pairs. The classical analog of the 
Gauss constraint can be satisfied but not in a unique way. Any pair that satisfy: 

apb -bpa = , (86) 
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will satisfy the Gauss constraint [175]. Therefore, it is useful to set a = p^ = in the 
sequel. There is still some residual gauge freedom but we shall not discuss it here. 
The co-triad can also be constructed as: 

sgnpc\pb\ ,„ , /I — 7 Pb -in /] — \Pb 



oj Tidx = ; — ^3 u-n: + v pc 1 — i'T2d'U — y\pc\-. — -TiSYn.0 dip . (ov) 

By comparison of (I85bl ) (or (|87] )) with (l84l ) one can see that the following identification may 
be made: 

Pb = VT{2M-T), Pc = ±T\ (88) 



Therefore, the degeneracy in (1871 ) at p^ = 7^ pc corresponds to the classical horizon 
whereas Pc = corresponds to the classical singularity. 

Next a basis is defined in the Hilbert space, these are denoted as 4= [|^, r) + |— ^u, r)] 
and are made up of eigenstates of the operators corresponding to pb and pc'. 

Pb 1^, r) = -jfi l/i, r) , Pc \p, t) = -jT l/i, r) . (89) 

The volume operator is also needed: 

V= fd^x Vldet^l = A7Ty^\\pb\ -^V = ^^T^AK\\Pb\ , (90) 

which is diagonal in the \^, r) basis and possesses eigenvalues: 

V^r = 27rj^/^\p\y^\ . (91) 

As well, the co-triad operator can be created. In the notation of [175], lOc := f" and 

ujb '■= sgnpby^lpcl- It is noted that the co-triad can be written in terms of the holonomy and 
volume: 

uJc = {27Tj)-'Tr{T^hii{h^\V}) , 

where Hr corresponds to the holonomy along an interval in the i?-direction. The operator 
version is given by: 

coc Im, r) = -i(2^7)-iTr (tsKr {h^\ y}) |^, r) = ^\fi\ (V|t + 1| - Vk - M) I/^, r) 

(92) 
Similarly, for uJb one can construct [175] 

^b\fJ', r) = ^sgn{p)y/\T\\fi, t) (93) 

As is usual, a general state can be expanded in terms of the above eigenstates: \iIj) = 

J2cfMT Ip-, t). 

Next the Hamiltonian constraint is constructed. One may pursue this in two ways. One 
way is to write the extrinsic curvature connection, K in (I71cl ). in terms of the modified 
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connection A. The second way is to regard K as the connection to be used. However, in 
this second case, the holonomies are to be constructed as functions of K, not A. It was 
shown in [175] that utilizing the second method results in a Hamiltonian constraint of the 
form: 



C-fffim — — K — =E E Cjik 

7' x/d^tCEj ' 



k o T?k 



YKb-^F 



(94) 



Here °F := dK + [K, K] and Q := dT = —sinOT^dO A dip (which can be calculated 
utilizing the triad associated with the standard two-sphere metric). Written in terms of the 
co-triad one has: 

e.^kT'E'^^E"'^ = - (47r7£(,))-^ e""^ °a;^ hi [4-', v] . (95) 



VWe) 

Here, h^ corresponds to the holonomy along an edge in the A;-direction (k = R, 6, (p) (the 
index k is summed over). £(^) = £(51) = £(^) = 6, where 6 is the length of the curve in the 
coordinate directions over which the holonomy is measured (superscripts o denote that the 
quantity is in the K connection representation). As well, the last part of the Hamiltonian 
constraint can be written in terms of the co-triad: 

where Arq = Ar^ = Ag^ = 5"^ and /i(jj) := hihjh~^h~^ . 

At this stage, all quantities are constructed out of holonomies (in °uj) and tetrads (in "uj 
and V) and therefore we can go to the quantum picture by simply replacing these with their 
operator analogues. In order to make the constraint Hermitian, the gravitational constraint 
is defined as C^^^ '■= \ C^ + C^ "^ . Without details, we quote the main results of this 
construction: The Hamiltonian constraint operator, acting on the states yields a difference 
equation for c^r- 

C'giv Cm, r =2(5 (y^\T+26\ + \/M j [C^ +2S,T + 2S-C^- 2S, r + 2s\ 

+ {^\r + 5\- yJ\T-5\) [(^ + 2,5)c, + 4.,. - (1 + 2-i^6^)^,c,,^ 
+ (// - 25)c^_45,r] + 2(5 [y[^^26\ + v1^) [c„ _ 2,, , _ 2, 

-C, + 25,r-2,]=^- (97) 

To analyze the behavior at the singularity, one starts at some positive value of r and utilizes 
Wh to evolve the c^ ,- to smaller values of r. The singularity resolution issue is insensitive 
to the choice of initial conditions (provided, of course, that they are not pathalogical). 

It turns out that the coefficients are always regular throughout the evolution for all values 
of r > as well as r < 0. From (1881 ) and (l89l ). r = corresponds to the classical singu- 
larity. Therefore, the evolution of the c^^,- coefficients is regular and may proceed beyond 
the classical singularity. In this theory of quantum geometry then, the quantum analogue of 
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T = is no longer a boundary of the space-time. In essence there is a smooth "bounce" 
which evolves to another large region of space-time. Expectation values of curvature encod- 
ing quantities are also finite. Interestingly, (though perhaps not surprising) mini-superspace 
reduced LQG cosmological models possess similar behavior where the big-bang is replaced 
by a smoothly evolving quantum bounce. (See [178] and references therein.) 

6 Black Holes in Astrophysics 

Although there exists no unquestionable proof that black holes exist in nature, there is 
mounting evidence suggesting that black holes are present in our universe. The evidence 
must be, by definition, indirect. We overview here some of the observational evidence for 
the existence of black holes as well as how researchers utilize observations to deduce the 
properties of these fascinating objects. 

Perhaps the first evidence that some extreme objects exist in the universe was the ob- 
servation of X-ray sources outside the solar system and of quasars in the 1960s [179]. 
Quasars are objects which possess luminosities on the order of 10^^ that of the sun. Matter 
accreting into a black hole could most easily explain such emissions of electromagnetic 
radiation. It is now also believed that black holes are related in some way to the observed 
gamma ray bursts. Since it is expected that almost all black holes have some amount of 
rotation, the Kerr solution provides a viable background in which to study these emissions. 
Astrophysical theory is suggestive that there is a limit on the angular momentum param- 
eter of — 0.998Af < a < 0.998Af [180] with high rotations amost certainly containing 
an event horizon and therefore unlikely to be an alternative to a black hole without event 
horizon [181]. 

One of the earliest (and brightest) X-ray sources detected was Cygnus X-1 [182], [183], 
which was noted to vary with time. A large number of X-ray sources have since been 
discovered, many of them associated with optically faint, distant stars. In such cases, it 
was not possible that the star itself could be emitting the X-rays. Instead, an argument was 
put forward that the X-rays originated from the accretion of the star's outer material onto a 
yet unseen companion object, likely a neutron star [184] [185]. Further, it was postulated 
that the slightest amount of angular momentum in the material (likely inevitable) would 
preclude anything resembling radial in-fall and the material would, rather, be forced into a 
disk around the compact companion (see figurefTTl). Viscous drag forces would then heat up 
the material to high enough temperatures to emit the radiation observed. The Uhuru satellite 
confirmed that the stars indeed must be orbiting some companion object, which must be 
very compact. If the mass of the companion is above the neutron star limit (approximately 
1.5-4 solar masses) then the likely alternative is a black hole. 

The evidence for a binary system comes from the periodicity of the visible star's spec- 
trum. One then needs to determine if the properties of the unseen companion allow it to be 
a neutron star or some more compact object. Consider a binary system as shown in figure 

As an approximation, Kepler's law can be utilized: 

2 _ Mi47r2ai 
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Figure 11 : A rotating black hole accreting matter from a nearby star Although the outskirts of the 
accretion disk is tilted with respect to the orbital plane, the inner regions are forced into alignment 
with the orbital plane of the black hole. 




Figure 12: An example of a binary system. 



where T is the orbital period and Mi and M2 represent the mass of the dark companion 
and observable star respectively. If the period, M2, and a2 are provided, the mass of the 
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the companion may be determined. In practice it is usually difficult to determine 02 as it 
requires knowledge of the inclination angle, 6i. 

For the case of Cygnus X-1, the optical member is a hot OB star (HDE 226868) which 
typically have very large masses (> 20 solar masses). With best estimates for the param- 
eters, the conclusion is that the dark companion in this binary system possesses a mass of 
the order 10 solar masses, much greater than the neutron star limit. It is therefore likely that 
the companion is a black hole. One of the sources of uncertainty is the value of the mass of 
HDE 226868. 

A class of X-ray objects are known as the X-ray transients. These objects emit X-rays 
periodically, followed by long periods of no emission. This allows one to study the optical 
companion in detail during the X-ray-quiet periods without the noise from the X-rays inter- 
fering with the observations. This is useful since detailed studies of the optical member of 
the binary will yield tighter constraints on the mass of the star and therefore on the mass of 
the companion object. In such a system, if the compact companion possessed a solid sur- 
face, it should be possible to see a characteristic emission of energy as the accreting material 
is brought to rest on the surface [186], [187], [188]. Observations of various sources seem 
to show no indication of such emission, and therefore the presence of an event horizon, as 
opposed to an object with a solid surface, is favored. 

As mentioned earlier, the discovery of quasars in the 1960s has led to speculation that 
some compact object must be responsible for the emission of so much energy. Over the 
years the evidence has become compelling that the gravitational sources are likely super- 
massive black holes at the center of galaxies. These objects are now generally referred to as 
active galactic nuclei (AGN). These sources typically emit 10^^ — lO^'* solar luminosities 
and have length scales on the order of less than a light-year and in many cases the scale may 
be measured in light-hours. These scales are based on the fact that appreciable changes in 
a system can not occur on time-scales shorter than it takes light to cross the system. 

The physics involved in AGNs is similar to the compact binary described above. 
Namely, nearby matter is accreted into a disk around the black hole and X-rays are emitted 
via a friction mechanism. A natural question arises in the case of AGNs: Could the grav- 
itational effects required to produce such X-ray emissions be due to the large number of 
stars and galactic matter near the core of the galaxy instead of a black hole? The constraints 
on the size, the lack of periodicities in the signals and the stability of the signals seem to 
favor a single central object (with masses of the order of lO^'' solar masses!) instead of a 
widespread, non-uniform source [189]. Also, there are "jets" of material present in may 
AGNs which remain aligned for time periods on the order of 10^ light-years (see figure [T3]). 
This indicates that a preferred axis must have been present in the system for at least that 
long, making a gravitationally bound compound object unlikely. 

Modeling these jets is an extremely difficult task involving general relativistic magne- 
tohydrodynamics. Large-scale computing must be employed in order to produce reliable 
results from the models. The jets can arise from a complex interplay between gas evapo- 
rating off of the accretion disk and magnetic fields present, known as the Blandford-Znajek 
process [190]. 
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Figure 13: Jets of ionized gas being ejected from the galactic core of M87. Permission to use 
image, made available through NASA and STScI, kindly provided by Tod R. Lauer (National Optical 
Astronomy Observatory), Sandra M. Faber (UC Observatories/ Lick Observatory), (figure quality 
reduced for arXiv file size) 



One of the major sources of the X-rays is the iron Ka line and there are several effects 
on this line. One effect is not strictly speaking gravitational in origin. It is the special 
relativistic doppler shift, where the line is blueshifted on one side of the accretion disk and 
redshifted on the other side. This will yield a doppler broadened line. Another effect is also 
present in the absence of gravity, this is the special relativistic beaming effect whereby the 
radiation intensity is amplified in the direction of particle motion compared to the intensity 
in the rest-frame. This effect is also differential in that the effect enhances the intensity 
on the approaching side of the disk. Another effect is strictly gravitational in origin. This 
is the gravitational redshift and time dilation. Unlike the doppler shift, this shift does not 
depend on what side of the disk (approaching or receding side) the material is residing. 
Also, the gravitational time dilation has the effect of reducing the overall flux since the 
emitter is "slowing down" compared to an observer at infinity. These effects act to skew 
the line profile, as is illustrated in[l4l The difference between the X-ray spectrum of matter 
accreting into a spinning versus non-spinning black hole is also displayed in this figure. This 
difference due to the spin arises from the fact that stable circular orbits in Schwarzschild 
geometry do not exist below r = 6Af , whereas for a rotating black hole this value is 
much smaller. Therefore, the gravitational redshift and time dilation can be much more 
pronounced in the case of a Kerr black hole as it is generally expected that the bulk of X- 
ray emission occurs in orbits at or above the stable orbit limit. With this assumption, data 
from the XMM-Newton satellite, analyzing the X-ray iron line of Seyfert galaxy MCG-6- 
30-15 constrains the Kerr spin parameter to be \a\ > 0.93 [191]. The rotational dragging of 
a Kerr black hole also has the effect of forcing the portion of the accretion disk that is close 
to the black hole to orbit in the equatorial plane (see figure fTTl). 
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Figure 14: An example of an X-ray spectrum of iron atoms in the vicinity of a non-spinning (left) 
and spinning (right) black hole. A spinning black hole drags the particles around it, via the Lense- 
Thirring effect, this allows for particles to orbit nearer to the black hole. Kind permission to display 
this figure, from the Harvard Chandra website, was granted by NASA/CXC/SAO. 



Interestingly, there is strong evidence that there is a super-massive black hole at the cen- 
ter of our own galaxy (Sgr A*). The "close" proximity of this black hole allows astronomers 
to directly measure its influence on its stellar neighbors [192] - [195]. The optical range of 
frequencies cannot penetrate the galactic center so studies of the galactic black hole are usu- 
ally performed utilizing the radio or the infra-red region. Knowing the orbital parameters 
between the black hole and its nearby stars, Kepler's law allows the determination of the 
approximate mass of the black hole. A mass of approximately 3.5 x 10^ solar masses is 
calculated for the mass of the galactic black hole [196], [197]. Observational data relating 
to these orbits, along with Keplerian fits, are displayed in figure [Tsl 

Other possible methods to detect black holes include: Gamma ray bursts, gravitational 
lensing (of background objects as well as of the orbital features of the accompanying binary 
star) [198], [199], and hopefully in the near future, gravitational wave signatures. Other 
proposals may be found in [200], [201] and [202]. 

We have only scratched the surface here regarding the observations and theoretical tech- 
niques used to study black holes in astrophysical contexts. There exist a number of excellent 
books and reviews on the subject. The interested reader is referred to the (much more thor- 
ough) review articles [203], [204] and [205] and the large number of references therein. 

7 Alternatives to Black Holes 

In this final section we will mention a few proposed alternatives to black holes along with 
possible measurements that may be performed in order to distinguish these objects from 
black holes. Some of these objects are black holes in the strict sense of the word. That is. 
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Figure 15: Observation of the orbits of stars near the galactic center This figure displays data 
points and best fit Keplerian orbits. Figure reproduced from [196] with kind permission from F 
Eisenhauer and The Astrophysical Journal, (figure quality reduced for arXive file size) 



they may contain a horizon but there is no singularity hiding behind it. It would be difficult, 
but not impossible, to distinguish some of these models from black holes [206]. We list 
some of the alternatives here with a brief description. 

Neutron stars with non-standard equation of state: Perhaps the greatest Achille's heel to 
the arguments in favor of black holes as the likely candidate in the binary systems discussed 
above is the uncertainty in the equation of state. This is the main reason for the large uncer- 
tainty in the neutron star limit. The regime of neutron star density is above what can reliably 
be studied in a lab and therefore any properties at these densities are not well constrained. 
Since pressure is also a source of gravity in general relativity, modification of the equation 
of state could increase the maximum mass that neutron stars may possess and therefore 
some large mass objects thought to be black holes could turn out to be neutron stars. A 
general form of the equation of state was studied in detail by Rhoades and Ruffini [207]. 
They made mild assumptions such as the speed of sound being bounded < c^ < 1, and 
that at lower densities (below some value po), it should produce equations of state thought 
to be well understood. In their study, they found a neutron star limit of approximately 3.2 
solar masses. Adding rotation to the picture yields [208]: 



M^ 



8.4 



10 



14 



-1/2 



M^ 



(99) 



It is therefore not inconceivable that the neutron star limit could be as high as 8-10 solar 
masses or higher. 

A possible scenario is the "Q-star", which allows for the possibility of nucleon confine- 
ment under extreme conditions. In these theories, it is expected that under certain condi- 
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tions, the equation of state differs strongly from standard ones, even at for relatively low 
densities (low values of po). Therefore, the assumptions used to derive the expression above 
are no longer valid. Such stars may possess masses on the order of 100 solar masses yet pos- 
sess radii which are approximately 1.4 times the corresponding Schwarzschild radius [209]. 
However, as mentioned previously, the absence of flare-ups due to material being 
brought to rest on a hard surface is in favor of a black hole instead of a neutron star or 
Q-star. As well, no reasonable value of po would allow a neutron star scenario for the 
super-massive galactic black holes thought to be responsible for AGNs. 

Repulsive interiors: This is not an alternative to a black hole as much as it is a possible al- 
ternative to the standard picture of an event horizon shielding a singularity. These scenarios 
basically stem from the fact that there is no reason to believe that the Schwarzschild solution 
is valid down to r = 0. In the T-domain (the r < 2M domain of the Schwarzschild solu- 
tion) it is possible, for example, to patch the Schwarzschild solution to a deSitter metric via 
a shell located at some space-Uke surface. The idea is to preserve the properties of the event 
horizon, which seems to fit observational data, but modify the interior. This idea seems 
to date back to Sakharov and Gliner who considered the possibility that, under extreme 
conditions, matter would possess an equation of state of the form p = —p [210], [211]. 
Explicit constructions of this model were performed in [212] and alternates of this model 
were also considered in [213]- [215]. A conformal diagram of Schwarzschild space-time 
with a deSitter interior is shown in figure [ 



r=oo 




deSitter universe 



Figure 1 6: A conformal diagram illustrating the space-time that results from patching the T-domain 
of the Schwarzschild space-time to a deSitter universe at a junction shell. Presumably a phase- 
transition of the collapsing matter occurs at the shell yielding the deSitter interior. T denotes the 
interior time coordinate whereas r denotes the exterior radial coordinate. 



Gravastars: A recent extension of the above idea is the gravitational vacuum star, or 
gravastar. The gravastar idea originated with P. Mazur and E. Mottola as an alternative 
to a black hole and possesses no event horizon [216] - [218]. In the gravastar picture, quan- 
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tum vacuum fluctuations are expected to play a non-trivial role in the collapse dynamics. A 
phase transition is believed to occur yielding a repulsive de Sitter core which aids in bal- 
ancing the collapsing object and thus preventing horizon (and singularity) formation. This 
transition is expected to occur very close to the limit 2m{r)/r = 1 so that, to an outside 
observer, it would be very difficult to distinguish the gravastar from a true black hole. The 
idea of a phase transition of the vacuum from a A ss state to a non-negligible A state is 
motivated from the behavior of Bose-Einstein condensates. The final gravastar configura- 
tion would also possess much less entropy than a black hole of similar size and therefore 
the problem of where the enormous black hole entropy comes from is alleviated. 

The original Mazur - Mottola model consisted of a deSitter interior separated from a 
Schwarzschild exterior via a finite shell with an equation of state satisfying p = +p (with 
thin shells on either side for patching purposes). It was later shown that, were a transition 
between a deSitter center and Schwarzschild exterior to be smooth and yield physically 
reasonable outer layers, anisotropic pressures must be present within the structure [219]. 
Models with continuous pressures satisfying various equations of state were explicitly con- 
structed in [220]. Examples of the pressure and density profiles for a sample gravastar 
(originally displayed in [220]) are displayed in figure [T7] Lobo and Arellano have studied 
several variants of gravastars or gravas tar-like objects in [221] [222]. 




3M 



Figure 17: The gravastar with energy density profile p = po exp[— (r/ro)"] and anisotropy A = 
a^ [p/ Po) A'/IS. Displayed are: radial (lower solid Une) and transversal (upper solid line) pressures, 
energy density (dashed line) and the compactness (dotted line). In this example the parameters 
are n = 3, total mass of configuration M = 1 and maximal compactness within the gravastar is 



/^max — 0.80. Notation is as follows: A 



2ra(r) p 



, with a and po constants, p is the 



Po 12 r Po ' 

"compactness" function 2m(r)/r. The transverse pressure, radial pressure and energy density are 
denoted as pt, Pr and p respectively. 



Note that by definition, the strong energy condition cannot be satisfied in any model 
with a deSitter region. A thorough discussion of gravastar energy conditions may be found 
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in [223]. 

In the case of anisotropic models, one way to distinguish their presence from either a 
black hole or a neutron star is via their surface redshift. It is known (see [224] and references 
therein) that stability in anisotropic spheres allows for a higher maximum redshift than for 
a stable perfect fluid sphere of similar mass, due to increased allowable compactness. Also, 
Chirenti and RezzoUa have discussed how to distinguish a gravastar from a black hole via 
quasi-normal mode analysis [225]. 

The list presented above is not exhaustive but should cover some of the most popular 
alternatives to black hoes. Another interesting alternative for example, put forward by 
Robertson and Leiter, is the magnetic eternally collapsing object (or MECO) (see references 
[226], [227], [228] and references therein.) 
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